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THE HOST GALAXY AND THE EXTENDED EMISSION-LINE REGION OF THE RADIO GALAXY 3C 79 1 
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ABSTRACT 

We present extensive ground-based spectroscopy and HST imaging of 3C 79, an FR II radio galaxy associated 
with a luminous extended emission-line region (EELR). Surface brightness modeling of an emission-line- 
free HST 7?-band image reveals that the host galaxy is a massive elliptical with a compact companion 0."8 
away and 4 magnitudes fainter. The host galaxy spectrum is best described by an intermediate- age (1.3 Gyr) 
stellar population (4% by mass), superimposed on a 10 Gyr old population and a power law (a\ = -1.8); the 
stellar populations are consistent with super- solar metallicities, with the best fit given by the 2.5 Z models. 
We derive a dynamical mass of 4 x 10 11 M within the effective radius from the velocity dispersion. The 
EELR spectra clearly indicate that the EELR is photoionized by the hidden central engine. Photoionization 
modeling shows evidence that the gas metallicity in both the EELR and the nuclear narrow-line region is 
mildly sub-solar (0.3-0.7 Z ) — significantly lower than the super-solar metallicities deduced from typical 
active galactic nuclei in the Sloan Digital Sky Survey. The more luminous filaments in the EELR exhibit 
a velocity field consistent with a common disk rotation. Fainter clouds, however, show high approaching 
velocities that are uncoupled with this apparent disk rotation. The striking similarities between this EELR and 
the EELRs around steep- spectrum radio-loud quasars provide further evidence for the orientation-dependent 
unification schemes. The metal-poor gas is almost certainly not native to the massive host galaxy. We suggest 
that the close companion galaxy could be the tidally stripped bulge of a late-type galaxy that is merging with 
the host galaxy. The interstellar medium of such a galaxy is probably the source for the low-metallicity gas in 
3C79. 

Subject headings: galaxies: active — galaxies: abundances — galaxies: individual (3C 79) — galaxies: inter- 
actions — galaxies: ISM 



1. INTRODUCTION 

Low-redshift quasars are often surrounded by massive ion- 
ized nebulae showing filamentary structures on scales of a few 
tens of kpcs. More specifically, about 40% of the quasars 
at z < 0.5 that are also steep- spectrum (i.e., usually FRII 
type) radio sources show such extended emission-line re- 
gions (EELRs) with an [O III] A5007 luminosity greater than 
4 x 10 41 erg s _1 . Here we will focus our discussion on these 
steep-radio- spectrum quasars, because (1) they are possibly 
the counterparts of FRII radio galaxies like 3C79, and (2) 
the presence of a powerful radio jet seems to be a necessary 
(though not sufficient) condition for producing a luminous 
EELR, as implied by the correlation between radio spectral 
index and extended optical emission dBoroson & Okel[l984l: 
IStockton & MacKentvl[l987h . 

An EELR typically has a mass of 10 9-10 M , displays glob- 
ally disordered kinematics, and shows a complex morphology 
that bears no obvious relationships ether with the host galaxy 

1 Based in part on observations obtained at the Gemini Observatory, 
which is operated by the Association of Universities for Research in As- 
tronomy, Inc., under a cooperative agreement with the NSF on behalf of the 
Gemini partnership: the National Science Foundation (United States), the 
Particle Physics and Astronomy Research Council (United Kingdom), the 
National Research Council (Canada), CONICYT (Chile), the Australian Re- 
search Council (Australia), CNPq (Brazil) and CONICET (Argentina). Gem- 
ini Program ID: GN-2006B-C-3. Some of the data presented herein were 
obtained at the W.M. Keck Observatory, which is operated as a scientific part- 
nership among the California Institute of Technology, the University of Cali- 
fornia and the National Aeronautics and Space Administration. The Observa- 
tory was made possible by the generous financial support of the W.M. Keck 
Foundation. Based also in part on observations made with the NASA/ESA 
Hubble Space Telescope, obtained from the Data Archive at the Space Tele- 
scope Science Institute, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under NASA contract NAS 5-26555. 



or wi t h the extended radio structur e (IStockton & M acKentv 
119871: iFu & Stockton! I2006L 12007 bt seeHtockton etaL|l2006 
for a review). In terms of quasar luminosity, host galaxy 
luminosity, and black hole mass, the quasars that show 
EELRs (hereafter the "EELR quasars") do not differ signif- 
icantly from non-EELR quasars. But the former are all low- 
metallicity quasars. The gas metallicity in their broad-line 
regions (BLRs) is significan tly lower (< 0.6 Zq) t han that of 
non-EELR quasars (> Z ) (IFu & Stocktonll2007al) . Although 
this gas-phase metallicity is unexpectedly low for galaxies 
with the masses found for these quasar hosts, it is consis- 
tent with that in their EELRs. Combining all of the pieces 
together, the most likely scenario for the origin of an EELR 
is that it comprises gas that is native to a gas-rich galaxy that 
merged with the quasar host and triggered the quasar activity, 
after which a large fraction of the gas was impulsively swept 
out by a large-solid-angle blast wave (i.e., a superwind), in a 
manner similar to that envis ioned for quasar-m ode feedback 
in the early universe (e.g., Di M atteo et aD2005l) . 

3C79(z = 0.256,l" = 4 kpc) 2 is a narrow-line ra dio galaxy 
havin g FRII radio je ts wit h a central compo nent (iHes et al.l 
fl995l: ISpangler et al.l fT984l: lHardcastle eTal][T997h . Hubble 
Space Telescope (HST) WFPC2 broad-band images show a 
complex optical morphology — a bright elliptical galaxy 
with an effective radius of 7.5 kpc, accompanied by multi- 
ple "ti dal arms" and two distinct "cores" othe r than the nu- 
cleus (iDunlopitan [200l; [de^ffetaDfllli). As we will 
show in § 13.11 only the galaxy and the "core" southwest 
of the nucleus are true continuum sources; the rest are all 
line-emitting regions. The nuclear spectra of 3C79 show 

2 Throughout we assume a flat cosmological model with Hq = 70 km s _1 
Mpc -1 , Q m = 0.3, and Q A = 0.7 
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Fig. 1 .— Hubble Space Telescope (HST) WFPC2 FR680N [O in] image of 
3C 79. Superimposed are contours from the continuum dominated WFPC2 
F675W image, at levels of 1.6 x l(T 17 x(l, 1.53, 2.34, 3.57) erg cnT 2 s" 1 
A -1 arcsec -2 above the background, showing the structures near the nucleus. 
Also overlaid are the field of the GMOS/IFU (white rectangular box) and the 
position of the LRIS slit (white dashed lines; PA = 142°). The extraction 
apertures for the Nl and SI regions are shown with grey boxes. The FRII 
radio jet directions are indicated by the black dashed lines. The angular size 
of the scale bar is 2"51. North is up and east is to the left for this and the 
following images. 



stellar absorption lines and a red continuum, indicating a n 
old stellar population dMillerl 1198 ll: iBoroson & 6kelll987l) . 
A line-emitting region in a curving filament extending 12" 
to the northwest of the galaxy was seen in an Ha image 
and was subsequently confirmed by a slit spectrum coverin g 
the [Oin] AA4959,5007 region (McCarthy et al. 1995, 1996). 
The spectrum also shows emission-line clouds between the 
filament and the nucleus, as well as clouds extending 5" to the 
south east. An [Q Hi] A5 QQ7 image taken in much better see- 
ing (IStockton et al.ll2006|) shows a very rich morphology and 
an emission-line luminosity comparable to that of the EELRs 
around steep- spectrum radio-loud quasars, establishing that 
3C 79 is associated with an EELR. 

In line with the orientation-based unificati on schemes o f 
FRII radio galaxies and radio-loud quasars (Barthel 1989), 
we believe that 3C 79 belongs to the same class of object as 
the EELR quasars, such as 3C 249.1 and 4C 37.43, but it is 
viewed at a different angle with respect to the molecular torus 
surrounding the central engine. If this is true, then the unique 
geometry of radio galaxies like 3C 79 could afford us a natu- 
ral coronograph blocking the blinding glare from the central 
engine. Such objects can be especially useful for studies of 
their host galaxies, which potentially preserve a large amount 
of information regarding the forces that have sculpted these 
spectacular emission-line regions. 

In this paper, we study the host galaxy and the EELR of 
3C 79 with extensive ground-based spectroscopy in combina- 
tion with a re-analysis of archival HST WFPC2 multi-band 
images. First, we briefly describe our observations and data 
reduction procedures in §12 In §[3] we study the morphology, 
stellar kinematics and stellar populations of the host galaxy. 
We then demonstrate the similarities between the 3C 79 EELR 
and the ones around quasars, in terms of gas kinematics, pres- 
sure, ionization mechanism, and metallicity in § [4] Finally, 
we discuss our main results in §[5] and close with a summary 
in§0 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. GMOS Integral Field Spectroscopy 

We observed 3C79 with the Integral Field Unit (IFU; 
lAllington-Smith et al. 2002) of the Gemini Multiobject Spec- 
trograph (GMOSi lHook et alJ2004h on the Gemini North tele- 
scope. The observations were performed on the night of 2006 
December 21 (UT). Three exposures of 2880 s were obtained 
using the half-field mode with the B600/G5303 grating at a 
central wavelength of 6242 A. With this setting, we had a 
field-of-view (FOV) of 3. // 5 x 5", a wavelength range of 4100 
to 6900 A, a dispersion of 0.46 A pixel -1 and an instrumental 
full width at half-maximum (FWHM) of 4.5 pixels (2.1 A). 
The host galaxy was placed at the lower-left corner of the IFU 
field so that a large fraction of the western part of the EELR 
could be covered (Figure [T]). Feige34 was observed for flux 
calibration. 

The data were reduced using the Gemini IRAF package 
(Version 1.8). The data reduction pipeline (GFREDUCE) 
consists of the following standard steps: bias subtraction, 
spectral extraction, flat-fielding, wavelength calibration, sky 
subtraction, and flux calibr ation. Cosmic ray r ejection was 
performed by L.A.Cosmic ( vanDokkum 2001) before run- 
ning the data through the reduction pipeline, with careful ad- 
justment of the parameters to avoid misidentification of real 
data. Spectra from different exposures were assembled and 
resampled to construct individual data cubes (x,y,A) with a 
pixel size of 0."05 (GFCUBE). To correct for the differential 
atmosphere refraction (DAR), we first binned the data cubes 
along the wavelength direction to increase the S/N ratio of the 
host galaxy; then the centroid of the galaxy was measured for 
each bin, and a polynomial was fit to the (x,y) coordinates, re- 
stricted by the fact that DAR only causes the centroids to vary 
along a straight line in (x,y) plane (at the parallactic angle); fi- 
nally we shifted the image slices at each wavelength using the 
solution from the polynomial fit. The three data cubes were 
then combined and binned to 0."2 pixels, the original spatial 
sampling of the IFU fiber-lenslet system, to form the final data 
cube. 

2.2. LRIS Long -Slit Spectroscopy 

We obtained optical long-slit spectroscopy of 3C 79 on the 
night of 2007 October 17 (UT) with the Low Resolution 
Imaging Spectrograph ( Ok e et all 1 19951) on the Keck I tele- 
scope. The Cassegrain Atmospheric Dispersion Compensator 
(Cass ADC 3 ) was used during the observation. We took two 
1200 s exposures through a l /; slit centered on the nucleus at 
a position angle of 142° (Fig.Q]). On the blue arm we used the 
600 groove mm -1 grism blazed at 4000 A, whilst on the red 
arm we used the 400 groove mm -1 grating blazed at 8500 A, 
offering, respectively, wavelength ranges of 3 100-5600 A and 
5100-8900 A and spectral resolutions of 4 and 7 A (FWHM). 
The spectra were taken at low airmass (~1.1), and the seeing 
was about 0."7 throughout the night. 

The spectra were reduced in the standard fashion us- 
ing IRAF tasks. The blue-arm spectra were wavelength- 
calibrated using observations of arc lamps, whilst for the 
red-arm spectra we used night-sky lines. Background sky 
was removed by fitting low-order cubic splines to the re- 
gions on each side of the object along the slit direction using 
BACKGROUND. One-dimensional spectra of 3C79 were 

3 http://www2.keck.hawaii.edu/inst/adc/docs/ 
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Fig. 2.— HST WFPC2 images of 3C 79: (left) the F675W image, (middle) the FR680N [O in] image rotated, shifted and scaled to match the F675W image, 
(right) the residual from subtracting the FR680N [O in] image off the F675W image. The contours, which are from the same images but smoothed by a Gaussian 
function with a FWHM of 3 pixels, show the central high surface brightness features. Contour levels are the same as in Fig.[T] 



first extracted with a wide (1 \") aperture, and they were flux- 
calibrated and corrected for atmospheric absorption using the 
spectrophotometric standard Wolf 1346 (taken 5 min after the 
nautical twilight earlier in the night). The calibrated spectra 
from the two arms agree perfectly in their overlapping region. 
Then, one-dimensional spectra of high S/N ratio were ob- 
tained by using a 2" extraction aperture. These spectra were 
multiplied by appropriate smooth curves so that their continua 
match those of the wide-aperture spectra. Finally, the blue and 
the red sides, after flux calibration and atmospheric absorption 
correction, were joined and binned to a common linear disper- 
sion (1.86 A per pixel) using S COMBINE. The line -of-sight 
Galactic reddening (A v = 0.421. ISchlegel et al.ll998h was cor - 
rected with a standard reddening curve (ICardelli et al.lll989h . 
We obtained an absolute flux calibration by s caling the spec- 
trum to thetfSrphotometry result from §[3J](^a = 1.8 x 10" 16 
erg cm" 2 s" 1 A" 1 between 6300 and 7300 A). 

2.3. HST WFPC2 Imaging 

To study the host galaxy morphology, we obtained 
WFPC2/WFC F675W (4x565 s) and FR680N (2x1300 s) 
images, and a WFPC2/PC F702W (2x 140 s) image of 3C 79, 
available from the archive of the HST. Th e two broad- 
band i mage s have been pr e viousl y presented by lDunlop et al.l 
(120031) and IdeKoff etal.1 (119961) . The Linear Ramp Filter 
(LRF; FR680N) image has a central wavelength of 6284 A 
and a bandpass FWHM of 82 A at the position of 3C 79, i.e., 
the filter is centered on the redshifted [O in] A5007 emission 
line. Images observed with LRFs are not flat-field calibrated 
in the calibration pipeline, so we flattened the LRF image with 
an F631N flat field reference image. 

Additional cosmic rays after the standard pipe line reduction 
were identified and replaced by L.A.Cosmic (Ivan Dok kum 
120011) . again with careful adjustment of the parameters to 
avoid misidentification of real data. 

3. THE HOST GALAXY 

3.1. Galaxy Morphology 

Though the WFC2 F675W image is dominated by con- 
tinuum radiation, Uf3 and [O ill] AA4959,5007 lines fall on 
the shortward wing of the filter profile, where the transmis- 
sions are about 42, 76, and 85%, respectively. To study the 
morphology of the host galaxy, one has to remove the line- 
emitting regions from the image. We thus rotated the WFC3 
LRF [Olll] image to the orientation of the F675W image, 



shifted and scaled the image to match the F675W image, 
and finally subtracted it off the wide-band image to obtain 
an emission-line free image (hereafter the "F675W*" image; 
Fig. [2]). A close companion 0."8 southwest of the major host 
galaxy is seen in the F675W* image, which we designate as 
3C 79A. The host galaxy itself is simply called 3C 79 in this 
section. 

3C79A is better sampled in the PCI F702W image, so 
we used this image to stud y its morphology . We used C. 
Y. Peng's GALFIT software dPeng et al.ll2002k to fit the two- 
dimensional galaxy profile, and we used an oversampled Tiny 
Tim 4 PSF. The result is shown in Fig. [3] As indicated by the 
best fit Sersic index of n = 1.23, the profile is close to an ex- 
ponential. As shown by the residuals, the exponential model 
clearly provides a better fit than the de Vaucouleurs model. 

We chose to use the F675W* image to study the morphol- 
ogy of 3C79. We converted the geometric parameters from 
the best fit exponential model of 3C79A (effective radius, 
axis ratio, and position angle [PA]) to the WFC image, and 
froze them; however, we allowed the position and magnitude 
of 3C79A vary in the modeling of 3C79. Figure [4] shows 
that 3C 79 can be described by a de Vaucouleurs model, as 
indicated by the high Sersic index. Table \T\ summarizes the 
morphology results. The PA of the minor axis of 3C 79 is 
aligned within 5° to the direction of the radio jets (as defined 
by the hot spots; see Fig.Q]). The companion galaxy 3C79A 
is located 0."84 (3.3 kpc) southwest of 3C79, and its PA is 
al igned within 7° to th e direction of the host galaxy. 

iDunlop et all (l2003h modeled the 3C 79 host galaxy with 
a de Vaucouleurs model (plus a concentric PSF to account 
for the narrow-line region rNLR ]). For compa rison, their re- 
sults are also shown in Table [T] IDunlop et al.l used the origi- 
nal F675W image, i.e., the emisson-line regions were not re- 
moved; thus they overestimated the brightness of the galaxy 
and had a poor estimate for the axis ratio and PA. We were 
able to reproduce their results using GALFIT when modeling 
the original F675W image. 

3.2. Stellar Kinematics 

We extracted the host spectrum from the final combined 
GMOS data cube using a 0."5 (2 kpc) radius aperture 5 . Al- 
though the GMOS spectrum shows a poorer S/N ratio in com- 

4 http://www.stsci.edu/software/tinytim/tinytim.html 

5 We chose such a small aperture because: (1) it yields the best S/N ratio 
compared to apertures of other sizes; and (2) it avoids hitting the edges of the 
FOV at short wavelengths 



4 



Fu and Stockton 



3C 79A 
PC1 F702W 


Pure Model 


de Vaucouleurs 

* 




Exponential 

• 

% 




Sersic n = 1 .23 

1" 


* 5(* 1 * 



Fig. 3. — Two-dimentional modeling of the close companion galaxy 
3C79A. The top left panel shows the original WFPC2/PC1 F702W image. 
Directly below are best fit de Vaucouleurs ("r 1 / 4 law"), exponential, and Ser- 
sic models convolved with the PSF. To remove the light from the host galaxy 
and the nuclear NLR, we fitted a de Vaucouleurs plus PSF, and the best fit 
models are shown together with the 3C79A models. To the right of each 
model are shown the smoothed residuals after subtracting the models. The 
residuals are displayed at a higher contrast. The top right panel shows the 
best fit Sersic model without convolution with the PSF, which gives the best 
realization of the gross morphology of the galaxy. The white circles in some 
panels are centered at the position of 3C 79 A and have a radius of 0"35. 



parison with the LRIS spectrum, it has a better spectral res- 
olution (2.1 A FWHM as opposed of 7 A), making it best- 
suited for measuring stellar kinematic s . We chose to use the 
stellar synthesis models of Vazdekis] (119991) to fit the data. 
Although these models have very limited spectra l coverage 
compared with those of iBruzual & Charlotl (l2003[) . they have 
a better spectral resolution (~ 1.8 A), which is similar to that 
of the GMOS spectrum, and they do cover the spectral re- 
gion where the spectrum show the highest S/N ratio (S/N ~ 
33 at rest-frame wavelengths of 4800 A < A < 5500 A). For 



the modeling we used M. Cappellari' s IDL program 6 imple- 
mentin g the pixel-fitting method of Cappellari & Emsellem 
( 2004). The program finds the best fit to the data by con- 
volving the model templates with a line-of- sight velocity dis- 
tribution (LOSVD). Since we used a Gaussian function to pa- 
rameterize the LOSVD, the program reports the mean stellar 
velocity (V*) and velocity dispersion (a*). In Figure[5]we plot 
our be st fit model against the data. Following Emselle m et al.l 
(2004), we included a multiplicative Legendre polynomial of 
degree 6 in the fit 7 to correct the model continuum shape. The 
kinematics results are insensitive to the assumed model age, 
although x 2 increases rapidly for models younger than 3 Gyr 
and the oldest population yielded the minimum \ 2 - Using 
a Monte-Carlo approach, we found that the 1 a uncertain- 
ties of V* and a* are both about 11 km s" 1 . These errors 
should be regarded as lower limits since they do not account 
for the effect of the template and continuum mismatch. As- 
suming a negligible amount of systematic rotation, the mea- 
sured velocity dispersion of ~236 km s -1 implies a virial mass 
of 5Ra 2 /G = 1 .3 x 10 11 M© within 2 kpc from the galaxy cen- 
ter. We estimated the luminosity-weighted velocity dispersion 
within ri/2 (7.2 kpc) to be <j* /2 = 218 ± 20 km s~\ with the 
aperture correction function for SAURQN elliptical ga laxies, 
(<t* /2 /<t?) = ( ri/2 / r )-o o66±o.035 dCappellari et al.ll2006h . The 

virial mass inside r x / 2 is then M\ l j 2 = 4.0+q I x 10 11 M©. 

3.3. Stellar Population 

To study the stellar population of the host galaxy, we have 
performed detailed continuum modeling to the 3C79 nu- 
clear spectrum obtained from LRIS. The host galaxy shows 
a strong UV excess (Fig. [6]) and a broad (FWHM ~ 6500 
km s" 1 ) Mgll A2798 line (Fig. |7]), but no broad Ua or H/3 
lines are seen. These charac teristics seem to resemble closely 
the ra dio galaxy Cygnus A (lAntonucci et al.|[l994l;l0gle et alJ 
1 19971) , implying a significant scattered quasar component. 
In addition, stellar absorption features are evident (Call K 
A3933, the G band A4300, and the Mgi b A5175, etc.), im- 
plying a dominating old stellar population (OSP). As a first 
guess, we tried to model the continuum with three compo- 
nents: (1) a nebular continuum (at T = 10,000, 15,000, and 
20,000 K), (2) a power law (F\ ex X ax ;— 0.4 < a x < -1.8 in 
accordance with the co mposite UV/optical quasar spectrum; 
IVanden Berk eTaII[200lh . and (3) an OSP (7 Gyr < Age < 10 
Gyr 8 , andZ = 0.2, 0.4, 1.0, and 2.5 Z©; Z© = 0.02). The same 
IDL program was used in the modeling as in the previous sec- 
tion, but this time we did not include polynomials to adjust the 
model continuum shape. Since the nebular continuum scales 
with the B aimer lines, we required that the H7 flux from the 
models to be within 5% of the measured flux. 

Figured shows the best fit three-component model (i.e., 
the model that produces the least x 2 ), which is a combination 
of a 9-Gyr-old 2.5 Z population, a nebular continuum at T 
= 15,000 K, and a X a power law with a = -1.4. This model 
fits both the continuum shape and most of the absorption fea- 
tures, however, it underpredicts the Call K line and the G 
band. This inadequacy can be removed once we introduce a 
fourth component — an intermediate- age stellar population 
(1 Gyr < Age < 3 Gyr, and Z = 0.2, 0.4, 1.0, and 2.5 Z©). 
Figure[6f? displays the best fit four-component model, which 

6 http://www.strw.leidenuniv.nl/~mcappell/idl/ 

7 Applying an additive Legendre polynomial of similar degrees led to con- 
sistent results. 

8 At z = 0.256 the age of the Universe is 10.4 Gyr. 
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Fig. 4. — Two-dimentional modeling of 3C79 host galaxy. The top left panel shows the WFPC2/WF emission-line free image. Directly below are best fit 
Sersic and de Vaucouleurs models convolved with the PSF. The superimposed contours in these panels are from the same images but smoothed by a Gaussian 
function with a FWHM of 3 pixels, and are at the same levels as those in Fig. ^ To the right of each model are shown the original residual from the model 
subtraction and a Gaussian smoothed version. The top middle and top right panels show respectively the data and the best fit model (without convolution with 
the PSF), both have been magnified by a factor of two and are displayed at a lower contrast. For the companion galaxy (3C 79A), we used an exponential model 
and have fixed its geometric parameters (r^, b/a and PA) to those obtained from modeling the WFPC2/PC1 image (see Fig.[3j. 
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FIG. 5. — Modeling the kinematics of the central part (r < 2 kpc) 
of 3C79 host galaxy. The GMOS/IFU rest-frame spectrum of 3C79 host 
galaxy is plotted in black, with spectral regions excluded from the modeling 
highlighted in blue (including nuclear emission lines and a CCD chip de- 
fect). Overplotted in red is the best fit instantaneous burst model of Vazdekis 
(1999), with the metallicity and age labeled in red. Also labeled are the mea- 
sured radial velocity of the stars relative to that of the NLR (z = 0.25632) and 
the stellar velocity dispersion. The residual after being offset by 8 x 10~ 17 
erg cm -2 s _1 A -1 is shown with diamonds, and the dotted line going through 
the residual indicates the zero level. 



successfully reproduces the Ca II K line and the G band ab- 



sorption. The global best fit to the data is given by stellar 
population models with a metallicity of 2.5 Z . In order to 
gauge how robust this metallicity result is, we compared the 
best fit four-component models of different metallicities. The 
best fit 1.0 Z model gives a generally good fit to both the 
absorption features and the overall continuum shape, except 
that it cannot reproduce the high S/N continuum shape around 
the Mg I b feature between 5000 and 5500 A (rest- frame), re- 
sulting in a 16% increase in the x 2 value relative to the 2.5 
Z model. Models with sub- solar metallicities all fail to give 
an adequate fit, especially in the continuum below 4300 A 
(rest- frame). We thus conclude that the stellar population in 
the host galaxy has a super-solar metallicity. A conservative 
range of metallcities that can give an adequate fit to the LRIS 
spectrum is between 1 and 2.5 Z . 

We have also attempted to use a young stellar population 
(< 0.1 Gyr) to replace the power law component, but the best 
fit model underpredicts the higher (> H7) B aimer lines if it 
correctly predicts the H7 flux (due to the B aimer absorption 
lines from the young population), and it gives a much poorer 
fit to the rest-frame UV continuum (2900A < A < 3800A), 
in particular, it overpredicts the Balmer limit break at 3650 A. 

If we adopt the average total stellar mass from the two best 
fit models shown in Fig. [6l then the stellar mass inside ri/ 2 
is Af* /2 = 2.6 x 10 11 M . The M* /2 /M^ r 2 ratio implies a dark 
matter fraction of ~ 34% within one effective radius, in agree- 
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Fig. 6. — Decomposing the nuclear continuum of the 3C79 host galaxy: (a) the best fit three-component model, and (b) the best fit four-component model. 
Spectral re gions that were ex cluded from the modeling (nuclear emission lines and imperfect sky subtraction regions) are highlighted in blue. Overplotted below 
(green) are Bruzual & Chariot ( 2003) stellar population models, X a power laws representing quasar scattered light, and reddened nebular continua (including 
Balmer lines > H7). The sum of the models (red) provides the best fit to the observed spectrum. Labeled are the metallicity, age and stellar mass for the stellar 
population models, the temperature of the nebular continuum models and the index of the power laws. The residual is shown with diamonds, and the dotted line 
going through the residual indicates the zero level. 
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Fig. 7. — Detection of a broad Mg II A2798 line in the nuclear spectrum 
of 3C 79. Dotted curve shows the LRIS spectrum, the solid curve is a best fit 
model of eight Gaussians plus a linear continuum, and the long dashed curve 
shows the broad (FWHM ~ 6500 km s _1 ) Gaussian profile of the Mg II line. 
Key narrow emission lines are labeled. 

ment with the dark matter fraction of nearby early-type galax- 
ies (|Cappellari et al. 2006). We caution that there are signifi- 
cant uncertainties in both the stellar mass derived from single- 
burst stellar synthesis models and the dynamical mass using 
the virial theorem. Hence, this straightforward calculation of 
dark matter content is subject to an error of at least a factor 
of two. We note that the absolute 7?-band magnitude 9 of the 
galaxy is about -23.2, whi ch is the same as that of the host 
galaxies of EELR quasars dFu & Stoc kton 2007a|). 

The compact companion galaxy 3C 79A is detected in the 
GMOS data cube. The comparison of the spectrum extracted 
from 3C 79 A and a background spectrum from a symmetric 
location with respect to the minor axis of the host galaxy 
shows a red continuum 10 , indicating an old stellar population 
for 3C 79A. This galaxy is reminiscent of the cl ose companion 
galax y \" away from the EELR quasar 3C48 dStockton et al.l 
120071) . which is also elongated towards the AGN nucleus. The 
S/N ratio of the spectrum does not allow a detailed modeling. 

4. THE EXTENDED EMISSION-LINE REGION 

The HST LRF [O III] image (F ig. ED) shows not only th e fil- 
ament seen in the Ha image of McCart hy et al.l ([1995), but 
also luminous emission-line "arms" at distances around 2. "5 
from the nucleus, and an "inner arc" \" to the west. To com- 
pare the EELR of 3C 79 with those around quasars, we deter- 
mined the total luminosity in the [O III] A5007 line within an 
annulus of inner radius 11.2 kpc and outer radius 44.9 kpc 11 
(Aoiii]). We found a luminosity of L[onrj = 7.8 x 10 42 erg s" 1 , 
close to that of the most luminous quasar EELR at z < 0.5, 
i.e., t hat of 4C 37.43, Lronri = 9.4 x 10 42 erg s" 1 (converted 
from Stockton & M acKenty ll987|) . As a reference, the least 
lu minous detected EELR i n the r adio-loud quasar subsample 
of IStockton & MacKentvl (119871) has a luminosity about 10 
times lower (PKS 2135-147, L [om] = 6.9 x 10 41 erg s" 1 ). 

9 Converted from the F675W ST system to R-bmd Vega system (-0.7 
mag), after applyi ng ^-correction (-0 .27) and passive evolution correction 
(+0.23), following lLabita et al.1 ( J200^ . 

10 The continuum sha pe did not c hange much after we corrected it for the 
"fixed-aperture effect" (Fu & Stockton 2007b). 

11 This is equivalent t o the 10 kpc/40 kpc annulus used by 
Stockton & MacKenty (1987) since they assumed an empty universe with 
Hq = 75 km s _1 Mpc -1 and qo = 0. 



4.1. Gas Kinematics 

Kinematics of the ionized gas can be measured from strong 
emission lines. Since a single 48 min GMOS exposure is 
enough to acquire a good S/N ratio in the [O III] A5007 line re- 
gion, we derived the velocity fields from the single data cube 
that has the smallest airmass (AM = 1.172) and the best see- 
ing. To study the velocity structure for the clouds in the inner 
arc that is only \" from the nuclear NLR, we subtracted the 
light from the NLR and the bright nearby cloud 0."5 east to 
the NLR by simultaneously fitting two Moffat PSF profiles at 
each wavelength in the wavelength range where both clouds 
are apparent 12 . After trials it was found that the best subtrac- 
tion results came from using the best fit Moffat profile from 
modeling the Feige 34 data cube over the same wavelength 
range, indicating both data cubes have the same resolution 
(FWHM ~ 0."6 at the [O III] region). We then fixed the Moffat 
parameters to those of the standard star and let the positions 
and magnitudes vary. Finally, we fixed the positions to the av- 
erage values across the wavelength range from the previous fit 
and performed the fit with only the magnitudes as free param- 
eters. The resulting PSF- subtracted data cube is presented as 
a set of velocity channel maps in Fig. [8] The velocity field is 
shown in Fig. [9] The velocity dispersion measurements have 
been corrected for the ao = 42 km s -1 instrumental resolution. 
Due to the imperfect PSF subtraction process, no useful kine- 
matics can be extracted from the PSF- subtracted data cube 
in the nuclear region bordered by the white no-data pixels in 
Fig. [9] Hence, we measured the velocity field in this region 
from the original data cube and inserted the results into the 
final velocity field. 

Most of the extended emission comes from two main fila- 
ments — the inner arc and the outer arm 2. "5 from the nu- 
cleus. Both filaments show similar velocity gradients, with 
the northern part approaching towards us and the southern 
part receding from us at velocities of ~ 200 km s" 1 . To- 
gether, these two filaments are consistent with a common 
disk rotation with a flat rotation curve. This sort of locally 
ordered kinematics of the brighter clouds is reminiscent of 
that of the three most luminous clouds southeast of 3C 249.1 
dFu & Stocktonll2006|) . However, the much fainter clouds out- 
side of the outer arc disrupt this simple picture by showing 
large approaching velocities. 

The velocity dispersion over most of the region is around 
50 to 120 km s" 1 . The region showing large a (150-200 km 
s -1 ) between the inner arc and the outer arm, although roughly 
coincident with the jet direction, in fact also has multiple ve- 
locity components that happen to contribute light to the same 
pixels. 

Overall, the velocity field of the 3C79 EELR is similar 
to those of the quasar EELRs, in the sense that all of them 
appear globally disordered but locally ordered, and they all 
show super-sonic velocity dispersions (sound speed c s ~ 17 
km s" 1 ) within the same ran ge between 50 and 130 km s -1 
(iFu & Stocktonll2006ll2007bk . 

4.2. Spectra of Emission-Line Clouds 

To increase the S/N ratio, we have combined GMOS spec- 
tra within various regions. The spectra of the clouds in the 
inner arc (r < \"2) are heavily contaminated by light from 
the nuclear NLR, making it difficult to measure their line in- 
tensities accurately. We therefore concentrated on the EELR 

12 The two clouds have essentially the same velocity, at z = 0.25632 ± 
0.00001. 
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Fig. 8. — GMOS/IFU [O in] A5007 radial velocity channel maps for the 3C 79 EELR. The central velocities, in km s -1 , relative to that of the quasar NLR, are 
indicated in the upper right corner of each panel. The red multiplication signs mark the nuclear narrow-line region, and the red plus signs indicate the brightest 
nearby emission-line cloud, both of which hav e bee n subtracted from the data cube to show the EELR more clearly. Solid ellipses show the extraction apertures 
of the five emission-line clouds discussed in § |4.2| in individual panels corresponding to their central velocities. In the last panel these apertures are shown all 
together as dashed ellipses, along with their names labeled nearby. 
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Fig. 9. — GMOS/IFU [O in] A5007 velocity field of the 3C 79 EELR. Inside the FOV there is a white no-data boundary near the nuclear region. The velocity 
field outside of the boundary is measured from a data cube where both the nuclear NLR and the bright emission-line cloud directly to the east were removed (see 
also Fig. [8j, while inside the boundary the data come from the original data cube. The fluxes inside the boundary were scaled down by a factor of 6 to show 
the high surface brightness peaks. Overlaid are contours from the HST WF3 LRF [Om] image smoothed by a Gaussian function with a FWHM of 1.5 pixels; 
contours are at levels of 1.8 x 10~ 15 x (1,2,4,20,40) erg cm -2 s _1 arcsec -2 . The three panels, from left to right, are line intensity, radial velocity (relative to that 
of the nuclear NLR) and velocity dispersion maps. Pixels are 0"2 squares. The red multiplication signs indicate the center of the host galaxy. 



clouds in the outer arm and beyond. In Fig. [8] we have shown 
the extraction apertures for the four EELR clouds. We have 
also extracted a NLR spectrum using a circular aperture with 
a radius of 0."3. In order to measure the key metallicity di- 
agnostic — the [Nil] A6584/Ha ratio — we extracted LRIS 
spectra from three regions along the slit: a 2" -wide aperture 
centered on the nucleus was used for the nuclear NLR, and the 
apertures used for the Nl and SI regions are shown in Fig.Q] 



Unlike those for the EELR clouds, NLR spectra should be 
corrected for stellar absorption features. We fitted the con- 
tinuum using the same stellar synthesis model templates as 
in § 13.31 and subtracted off the best fit model to obtain the 
clean emission-line spectrum of the NLR. Note that the line 
ratios for the NLR should be reliable even if the continuum 
subtraction were just barely correct, given the large equiva- 
lent widths of the emission lines that are affected by strong 
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Fig. 10. — Log-scaled density distributions of SDSS Seyferts (blue contours), LINERs (red contours), and star-forming galaxies (including starburst-AGN 
composite galaxies; the background image) in line-ratio diagrams: (a) [O n] AA3726, 3729/H/? vs. [O in] A5007/H/3, and (b) [O n] AA3726, 3729/[Ne in] A3869 
vs. [O Hi] A 5007/H/3. Measurements fr om the EELRs of 3C79, 3C 249.1, and 4C 37.43 are shown as red open circles, violet solid circles and blue squares, 
respectively ( Fu & Stockton 2006, 2007b). Both diagrams clearly show that the EELRs are photoionized by an AGN-type ionizing spectrum. 



absorption features. Intrinsic reddening due to dust in front 
of or associated with the cloud was determined from the mea- 
sured H7/H/3 ratio for GMOS spectra and from Ha/H(3 ra- 
tio for LRIS spectra, assuming theoretical ratios of 0.468 for 
H7/H/3 (the val ue for case B reco mbination at T = 10 4 K and 
N e = 10 2 cm" 3 ; IOsterbrocklfT989h and 3.1 for Ha/Hp. For 
clouds a and b, where the S/N ratio of the H7 line is poor, 
we assumed zero intrinsic reddening. For clouds Nl and SI, 
where the measured Ha/H[3 is slightly lower than the theo- 
retical value, we also assumed zero intrinsic reddening. Both 
the intrinsic reddening and Galactic redde ning were correcte d 
using a standard Galactic reddening law (ICardelli et al.ll 989). 

The reddening-free emission-line fluxes and 3 a upper lim- 
its are tabulated in Table [2] The quoted 1 a errors were de- 
rived from the covariance matrices associated with the Gaus- 
sian model fits, with the noise level estimated from line-free 
regions on either side of an emission line. 

There is enough overlapping area to show that the data from 
the two instruments agree fairly well: (1) The low intrinsic 
reddening for cloud c is consistent with the measured Ha/H[3 
ratio of 3.09 in region Nl; (2) The GMOS NLR spectrum 
suggests an intrinsic reddening of Ay = 0.28, whilst the LRIS 
NLR spectrum shows Ay = 0.36; (3) The individual line ratios 
relative to H(3 of regions c and NLR (GMOS) are very much 
consistent with those measured from the LRIS Nl and NLR 
spectra, respectively. 

4.2.1. Electron Density and Temperature 

The [Oil] luminosity-weighted average electron density 
(N e ) and electron temperature (T e ) can be measured from 
the intensity ratio of the [On] AA3726,3729 doubl et and the 
rO m l (A4959 + A5007)/A4363 ratio, respectively (lOsterbrockl 
119891) . Our GMOS data cube covers all of these five lines, 
allowing us to measure the density and temperature in indi- 
vidual emission-line clouds. The [S II] AA6717,6731 doublet 
is covered by the LRIS spectra, however, they fall into the 
region where strong OH night sky lines dominate and severe 
CCD fringing occurs, preventing reliable measurements. We 
thus only discuss results from the [O II] doublet. 

The [O II] doublet is barely resolved in our spectra. How- 
ever, for the NLR and the EELR clouds c & d, the high S/N 
ratio profile of the doublet allows a meaningful decomposi- 
tion once we constrain their expected wavelengths to be at the 



same redshift as that of the nearby [Nem] A3869 line. Un- 
certainties of the line fluxes from this fit tend to be low as 
the wavelengths are no longer free parameters. Lifting this 
wavelength constraint gives a more conservative error esti- 
mate. Here we have adopted the mean values of the error es- 
timates from the two fits — one with and one without wave- 
length constraints. The weak [O in] A4363 line is detected 
in all the clouds except b. We then derived T e and N e con- 
sistently using the IRAF routine TEMDEN for the three re- 
gions where both [O II] and [O in] ratios are available. Ta- 
ble [3] summarizes the results. For the clouds in the EELR, the 
temperatures are about 12,000 to 14,000 K, and the densities 
are about 100-200 cm -3 , similar to those of quasar EELRs 
that were measured using the same technique dStockton et al.l 
12001 iFu & Stocktodl200ai2007bh . 

4.2.2. Ionization Mechanisms 

The emission-line ratios of the 3C 79 clouds are very sim- 
ilar to those of the quasar EELRs. Previously we have con- 
cluded that the spectra of quasar EELRs are most consistent 
with being photoionized by quasar continuum, and are incon- 
sistent with shock or the self- ionizing "shock + precursor" 
models dFu & Stock ton 2007b). Now we re-emphasize this 
point by comparing strong line ratios of the EELRs of the 
radio galaxy 3C79 and the quasars 3C 249.1 and 4C 37.43 
with emission-line galaxies at a similar redshift range. The 
Sloan Digital Sky Survey Data Release 4 (SDSS DR4) pro- 
vides high-quality optical spectra for a large number of galax- 
ies. The emission-line fluxes after correcting for stellar ab- 
sorption and foreground extinction are publicly available 13 . 
We first created a subsample of SDSS emission-line galaxies 
by including only those with S/N > 3 in the strong emission- 
lines [On] A3727 (=A3726 + A3729), H/5, [Om] A5007, Ha, 
[Nil] A6584 and [Sn] AA6717,6731. Then we divided the 
sample into star-for ming galaxies and AGN u sing the empir- 
ical dividing line of iKauffmann et al] (l2003|) . and corrected 
for the intrinsic reddening using the B aimer decrement an d 
the standard Galactic reddening curve (ICardelli etal.lH989h . 
The theoretical values for the Ha/H(3 ratio are 2.85 for star- 
forming galaxies and 3.1 for AGN. Finally, we classified 

13 h ttp://www.mpa-garching.mpg.de/SDSS/DR4/; and see Tremonti et al. 
( 2004) for a description of the data. 
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the galaxies into Seyferts, LINERs, star-forming galaxies 
and AGN/starf orming composites, following the classifica- 
tion schemes of Ke wley etail (l2006|) . 

Figure[l0] compares the locations of our EELRs and the dis- 
tributions of the SDSS galaxies in two line-ratio diagnostic 
diagrams. As can be seen, the EELRs share the same ter- 
ritory as the Seyfert galaxies, but are clearly distinguishable 
from star-forming and composite galaxies, indicating that the 
EELRs are photoionized by an AGN-type ionizing spectrum. 

The line-ratios from the 3C79 EELR match those from 
3C 249.1 and 4C 37.43 very well, thus pure shock and "shock 
+ precursor" models can be ruled out, since it has been shown 
(using the same line-ratios) t hat these models do no t fit these 
quasar EELRs (see Fig. 5 of iFu & Stockto n 2007b for an ex- 
ample). 

4.2.3. Gas Metallicity 

With the knowledge that the EELR of 3C 79 is photoionized 
by the hidden quasar, we can now estimate the metallicity of 
the gas using photoionization models. In the following we 
will compare our data with the du sty radiation-pressur e dom- 
inated photoionization models of IGroves et al.l (l2004|) with a 
density of 1000 cm -3 . Since the diagnostic line ratios we will 
use are almost entirely dependent on the abundances, models 
adopting a different density will not change our results. 

As the metallicity of a photoionized nebula increases, the 
temperature of the nebula decreases (as a result of the in- 
creased metal cooling), and the relative strengths between hy- 
drogen lines and lines from other heavier elements also de- 
creases (due to the increased abundance of the heavy elements 
relative to hydrogen). The latter effect is especially strong for 
nitrogen lines since the element is dominated by "secondary" 
production 14 at Z > 0.2 Z . Hence, N abundance increases 
at a much faster rate with metallicity, i.e., N/H ex Z 2 . In Fig- 
urefTTb we use the temperature sensitive [O in] A4363 / A5007 
ratio to demonstrate the dependance of nebula temperature 
on metallicity. The stratification among the grids of dif- 
ferent metallicities is further improved by the fact that he- 
lium increases lin early with metallicity relative to hydrogen 
dPagel et alJTl992h , In Fig. Gj}? and Fig. Q21 we show the 
strong dependance of [N I] A5 199 /Up and [N n] A6584/Ha on 
metallicity. At each metallicity, we show models spanning 
a range of ionization parameters (-2.3 < log(£/) < 0) and 
four power-law indices from a = —1.2 to -2.0, representing 
the quasar ionizing continuum (F u ocis a ). We have converted 
the modeled total metallicities to gas-phase metallicities using 
the solar abundances defined by lAnders & Grevessd J 19 89) 
(12+log(O/H) = 8.93 and Z = 0.02), in order to be con- 
sistent with both previous studies on the metallicity of q uasa r 
BLRs and metallicity re sults from stel l ar pop ulations (§ 13. 3b . 
One solar metallicity in IGroves etall (|2004|) corresponds to 
1/3 Z in our figures, since approximately half of the metals 
are assumed to be depleted onto dust in the models, and their 
assumed solar metallicity i s about 0.2 dex 15 lower than that of 
lAnders & Grevessd (119891) . 

As can be seen in Figs. [TT] & [121 the line-ratios of the 
EELR as well as the NLR of 3C 79, like those quasar EELRs 
and their BLRs, are most consistent with a low-metallicity 

14 where nitrogen is synthesized from existing carbon and oxygen 
(Pagel & Edmunds 198J). 

15 The C, N, O, and Fe abundances in the solar abundan ce set assumed by 
Groves et al. ( 2004) are all ~0.2 dex lower than those in Anders & Grevessd 
(1989). Hence, the mass fractio n of metals (i.e. Z) is ab out 1.5 times lower 
— Z© = 0.013 — than that of Anders & Grevesii tl989h . 



at about 1/3 to 2/3 solar. Fig. [121 illustrates the peculiarity 
of 3C79 more clearly by comparing it directly with SDSS 
AGN on the same line-ratio diagnostic. The nuclear NLR 
along with the two other regions where we have [N II] mea- 
surements all lie on the lower metallicity side of the main 
Seyfert branch, overlapping with quasar EELRs. This result 
supports the orientation-based unifica tion schemes of FR II 
radio galaxies and radio-loud quasars (Barthel 1989). On the 
other hand, 3C 79 offers yet another example strengthening 
(1) the correlation between the presence of EELRs and the 
metallicity of the nuclear gas (BLRs and/or NLRs), and (2) 
the similar metallicity in EELRs and the gas in the nuclear re- 
gion. Both correlations were initially discovered only a mong 
steep- spectrum radio-loud quasars (IFu & S tockton 2007 a|). 

5. DISCUSSION 

5.1. The Origin of the Low -Metallicity Gas 

The great majority of quas ars show super- sola r metallicities 
in their nuclear regions (see Ha mann et al .112007]. for a recent 
review ); yet, as we have previously shown dFu & Stockton! 
2007a), quasars with luminous EELRs are drawn exclusively 
from the subset of steep-radio- spectrum quasars that have 
BLRs with sub-solar metallicities. The fact that the quasar 
host galaxies are almost always very massive and thus ex- 
pected to have high metallicities suggests an external source 
for the low-metallicity gas. Furthermore, the apparent link 
between the metallicity of the gas in the BLRs and the metal- 
licity of the gas in the EELRs, which have minimum masses 
of 10 9-10 M , means that the external source of the gas must 
itself have been moderately massive. 

In 3C79, we have an FRII radio galaxy that shows the 
same pattern. Although we of course cannot obtain a broad- 
line-region metallicity in this case, we have used the unre- 
solved nuclear narrow-line region metallicity as a surrogate. 
Aside from providing a consistency check on unified models 
for FRII radio galaxies and quasars, the main question of in- 
terest is whether the clearer view of the inner region of the 
host galaxy gives us any insight into the origin of the low- 
metallicity gas and the mechanism that produced the EELR. 

The host galaxy of 3C 79 appears to be a fairly normal el- 
liptical morphologically. The most intriguing feature in the 
inner region is the extremely compact galaxy 0."8 from the 
center of the host galaxy. The nearly exponential profile of 
this galaxy suggests that it could be a tidally stripped "pseu- 
dobulge" of a late-type galaxy (e.g.. lKormendy & Kennicutt 
2004); the interstellar medium of such a galaxy would be 
quite plausible as a source for the l ow-metallicity gas in 
3C79. We have su ggested elsewhere (Stock ton et al.l 120071: 
IFu & Stoc kton 2007b) that luminous EELRs may result from 
nearly spherical blast waves connected with the initiation of 
FRII radio jets. The connection with low-metallicity gas 
is less clear, although we can speculate that it may have 
something to do with the lower radiative coupling such gas 
would have to the quasar radiation field, allowing more effi- 
cient accretion. It is interesting to note that related consid- 
erations have been invoked recently to explain the link be- 
tween low metallicit y and long-duration gamma-ray bursts 
dFruchter et al.ll2006|) on the scale of massive stars, although 
the physical mechanisms involved are certainly different. 

5.2. Low -Metallicity Radio Galaxies in SDSS 

Arm ed with the latest photoionization models. [Groves et aTl 
(120061) identified -40 (out of - 23000) candidates of low- 
metallicity Seyfert2s in SDSS. A caveat to note is that 
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Fig. 11.— Metallicity-sensitive line-ratio diagrams: (a) He II A4686/H/3 vs. [O ill] A4363/[0 ill] A5007, and (b) [Ni] A5199/H/3 vs. [O III] A5007/H/3 
Measureme nts from the EELRs of 3C 79, 3C 249.1, and 4C 37.43 are shown as red open circles, violet s olid circles and blue squares, respectively (Fu & Stockton 
2006, 2007b). Besides 3C79, accurate [Ni] flux measurements are only available for 4C 37.43 El (Stockton et al. 2002). Photoionization model grids from 
Groves et al. ( 2004) are shown for five metallicities from 0.1 to 1.3 Z©. For each metallicity, model predications are given for a range of ionization parameters 
(-2.3 < log(t/) < 0) and four power law indices representing the quasar ionizing continuum (F u oc v au ; ot v = -1.2,-1.4,-1.7, and -2.0). The grids for the 
highest metallicity are highlighted to show the inconsistency between our data and predictions assuming super- Solar metallicity gas. Both diagrams show that 
the EELR of 3C 79 consists of sub-solar metallicity gas, matching th e gas metallicity of t he EELRs of 3C 249.1 and 4C 37.43 (Fu & Stockton 2007b), as well as 
the metallicity of the BLRs of radio-loud quasars showing EELRs ( Fu & Stockton 2007a). 
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Fig. 12. — The nuclear NLR and EELR of 3C79, along with quasar 
EELRs, are much metal-poorer than in typical AGN, as seen in the metal- 
licity diagnostic diagram of [Nil] A6584/Ha vs. [O III] A5007/H/3. Keys 
are the same as in Fig. [TT] The background image shows the density dis- 
tribution of the SDSS AGN sample in log scale. Objects above the upper 
green curve are objects dominated by AGN (LINERs ar e concentrated in th e 
lower denser branch and Seyferts in the upper branch; Kewlev et al. 2006), 
below the lower curve are starforming galaxies (Kauffmann et al. 2003), and 
AGN/starforming composite galaxies are in between. Most Seyferts can be 
fit quite well with the super-solar metallicity models. The top region bor- 
dered by the black dashed line and the upper green curve is where the low- 
metallicity Seyferts are located ( Groves et al. 2006). 

the authors restricted the sample to galaxies with stellar 
masses lower than 10 10 M while selecting their candi- 
dates, as guided by the mass-metallicity correlation (e.g. 
iTremonti et al] l2004h . Although the EELR quasars are ex- 
cluded from the SDSS emission-line galaxy sample, 3C79 
represents a population of low-metallicity AGN that were 
missed by iGroves et all These AGN are hosted by massive 
evolved galaxies with ~ 10 12 M of stars and harbor ~ 10 9 
Mq black holes at their hearts, and they are likely to be more 
powerful than the Seyfert2s. 

By cross-correlat ing the TEXAS radio galaxies 
(iDouglas et al.1 Il996h with SDSS DR4 AGN, we have 
identified a sample of low-metallicity radio galaxies. Like 



that of 3C79, the spectra of their host galaxies show red 
slopes and absorption features that are indicative of an 
old stellar population. Follow-up high-resolution imaging 
and spatially resolved spectroscopy of these galaxies and a 
comparable control sample are needed to finally nail down 
the source of the low-metallicity gas and the triggering 
mechanism for quasar superwinds. In addition, detecting 
EELRs around these objects would provide a new test for the 
unification schemes. 

6. SUMMARY 

Based on extensive ground-based spectroscopy and archival 
HST WFPC2 images of the radio galaxy 3C 79, we have con- 
ducted a detailed analysis of its host galaxy and the EELR. 
The host galaxy of 3C 79 is a massive elliptical with M R = 
-23.2. The UV/optical spectral energy distribution of the host 
galaxy is best described by a combination of an intermediate- 
age stellar population (1.3 Gyr), an old stellar population (10 
Gyr), a power law, and a nebular thermal continuum. Both 
stellar populations are metal rich (2.5 Z ). This best fit model 
indicates a total stellar mass within one effective radius of ap- 
proximately 2.5 x 10 11 M , consistent with the virial mass 
derived from stellar kinematics (~ 4 x 10 11 M ). 

The EELR of 3C 79 is remarkably similar to the most lumi- 
nous quasar EELRs. The velocity field, although only avail- 
able in a small area, is locally ordered but globally disor- 
dered. The EELR is almost certainly photoionized by the hid- 
den quasar, and it shows densities of ~ 100 cm -3 in the [O II] 
emitting regions and temperatures around 13,000 K. Most in- 
terestingly, the metallicity of the gas in both the EELR and 
the NLR is about 1/3-2/3 solar, matching perfectly the metal- 
licity in both the EELRs and the nuclear BLRs of the EELR 
quasars. 

There is a very compact close companion galaxy 3.3 kpc 
from and 4 magnitudes fainter than the host galaxy. This 
companion galaxy shows an exponential profile, and is pre- 
sumably a tidally stripped "pseudobulge" of a late-type galaxy 
that might be the source for the low-metallicity gas as well as 
the starforming regions in 3C 79. 

The exact correspondence between this EELR and the 
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EELRs around quasars joins the already overwhelming evi- 
dence in support of the unification schemes for FRII radio 
galaxies and radio-loud quasars. The definitive trait of sub- 
solar metallicity in the NLR of 3C 79 also provides an effi- 
cient tool for selecting FRII radio galaxies that likely host 
luminous EELRs. 
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TABLE 1 

Morphology of 3C 79 Host Galaxy and Its Close Companion 



Galaxy 


Model 


F675W a 


F702W fl 


r 1/2 /kpc 


n b 


b/a 


P.A. 


3C79A 


Sersic 




22.19 


0.22 


1.23 


0.49 


18.2 




Exp 


22.24 


22.23 


0.22 


(1) 


0.48 


18.0 


3C79 


Sersic 


18.06 




10.1 


5.21 


0.70 


6.2 




de Vauc 


18.24 




7.2 


(4) 


0.71 


5.7 


Dunlop c 


de Vauc 


17.46 




7.5 


(4) 


1.0 


13 



a ST magnitudes corrected for Galactic extinction. 15 Sersic indexes. Converted from 
Dunlop et al. 1 2003) Table 3 to our cosmology. 



TABLE 2 

Line Ratios of 3C 79 Line-Emitting Clouds Relative to H/3 



Region 


[Ne V] A3426 


[O II] A3727 


[Ne III] A3869 


[O III] A4363 


HeIIA4686 


Hf3 


[OIII] A5007 


[NI] A5199 


Ha 


[Nil] A6584 


Gemini GMOS/IFU 


a 


0.83±0.19 


2.45 ±0.16 


0.65 ±0.07 


0.15 ±0.05 


0.33 ±0.04 


1.00 ±0.05 


9.32 ±0.03 


<0.18 






b 


1.54 ±0.43 


3.25 ±0.35 


1.18±0.19 


<0.49 


<0.33 


1.00±0.11 


8.91 ±0.08 


<0.47 






c 


0.53 ±0.10 


2.50±0.16 


1.01 ±0.05 


0.13 ±0.03 


0.26 ±0.02 


1.00 ±0.02 


10.39 ±0.02 


0.06 ±0.02 






d 


1.38±0.15 


2.62 ±0.09 


1.19 ±0.05 


0.19 ±0.02 


0.29 ±0.02 


1.00 ±0.02 


11.52±0.01 


0.05 ±0.02 






NLR 


1.057 ±0.016 


1.590±0.015 


0.944 ±0.004 


0.176 ±0.003 


0.302 ±0.003 


1.000 ±0.004 


10.923 ±0.002 


0.041 ±0.004 
















Keck LRIS/Lon^ 


*-Slit 










Nl 


0.60 ±0.02 


2.73 ±0.02 


1.09 ±0.02 


0.15 ±0.03 


0.28 ±0.02 


1.00 ±0.02 


10.63 ±0.02 


<0.09 


3.09 ±0.03 


1.24 ±0.03 


SI 


0.58 ±0.04 


2.92 ±0.05 


1.31 ±0.04 


0.27 ±0.07 


<0.15 


1.00 ±0.04 


11.83 ±0.03 


<0.21 


2.93 ±0.07 


0.91 ±0.06 


NLR 


1.093 ±0.003 


1.591 ±0.015 


0.921 ±0.009 


0.185 ±0.008 


0.297 ±0.007 


1.000 ±0.008 


10.602 ±0.008 


0.038 ±0.017 


3. 100 ±0.008 


1.040 ±0.007 



TABLE 3 

Properties of 3C 79 Line-Emitting Clouds 



Region 


V 


cr 




H/3 x 10 17 


[On] 


[Oin] 


N e 


T e 




(km s" 1 ) 


(km s _1 ) 


(mag) 


(erg cm -2 s _1 ) 


3726/3729 


(4959+5007)/4363 


(cm" 3 ) 


(K) 


a 


-145 


102 


(0.000) 


4.7 ±0.2 




85 ±29 






b 


-138 


91 


(0.000) 


1.4 ±0.2 




>23 






c 


-74 


109 


0.016 


21.0±0.4 


0.76 ±0.08 


108 ±22 


100 ±100 


12500^0 


d 


+131 


132 


0.392 


23.5 ±0.4 


0.80 ±0.08 


82 ±10 


150 ±100 


13900+800 


NLR 





136 


0.284 


99.6 ±0.4 


0.70 ±0.07 


83.0±1.4 


30 ±80 


13880 ±100 



a Intrinsic reddening. The values in parentheses are not directly measured from the H7/H/3 ratio, since the H7 lines in these spectra are not 
well-detected; and we assumed zero reddening for these clouds. See § 14.2. ll for details. 



